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ABSTRACT 

Recent work suggests that rich clusters of galaxies commonly have large populations 
of dwarf (ie. low luminosity) members, that is their luminosity function (LF) turns up to 
a steep slope at the faint end. This population, or more particularly the relative numbers 
of dwarfs to giants, appears to be very similar for clusters of similar morphology, but 
may vary between cluster types. We have previously suggested that dwarfs may be 
more common in less compact, spiral rich clusters. Similarly we have found evidence for 
population gradients across clusters, in that the dwarf population appears more spatially 
extended. In the present paper we summarise the current evidence and propose, in 
analogy to the well-known morphology - density relation, that what we are seeing 
is a dwarf population - density relation; dwarfs are more common in lower density 
environments. Finally we discuss recent semi-analytic models of galaxy formation in 
the hierarchical clustering picture, which may give clues as to the origin of our proposed 
relation. 



Subject headings: galaxies: clusters: general — galaxies: luminosity function, mass 
function — galaxies: photometry 
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1. Introduction 

Much recent work has been devoted to the 
question of the galaxy luminosity function (LF) 
within rich clusters, particularly with regard to 
the faint end which has become accessible to de- 
tailed study through various technical and ob- 
servational improvements (see e.g., Driver et al. 
1994; Biviano et al. 1995; Bernstein et al. 1995; 
Mohr et al. 1996; Wilson et al. 1997; Smith, 
Driver & Phillipps 1997 = Paper I; Trentham 
1997a,b). 

For the most part these studies concur that 
the LF becomes steep (Schechter (1976) slope 
a < —1.5) faintwards of about Mb = —17.5 or 
Mr ~ —19 (for H = 50 km s" 1 Mpc" 1 ), and 
Paper I suggested that such a dwarf rich popula- 
tion might be ubiquitous. In a subsequent paper 
(Driver, Couch k Phillipps 1997a = Paper III) 
we have examined the luminosity distribution in 
and across a variety of clusters, examining the 
possible dependence of the dwarf population (in 
particular the ratio of dwarfs to giants) on clus- 
ter type and position within the cluster. In the 
present paper we summarise the evidence to date 
for the (dis) similarity of the dwarf population in 
different environments. 

2. The Dwarf Luminosity Function 
2.1. Dwarfs in Rich Clusters 

Several papers (e.g., Driver et al. 1994; Pa- 
per I; Wilson et al. 1997) have recently demon- 
strated remarkably similar dwarf populations in 
a number of morphologically similar dense rich 
clusters like (and including) Coma. This simi- 
larity appears not only in the faint end slope of 
the LF, around a = —1.8, but also in the point 
at which the steep slope cuts in, Mr ~ —19 (ie. 
about M* + 3.5). The latter implies equal ratios 
of dwarf to giant galaxy numbers in the different 
clusters. 

However, there clearly do exist differences be- 
tween some clusters. For example, several of the 



clusters in the Paper III sample do not show a 
conspicuous turn up at the faint end. Either 
these clusters contain completely different types 
of dwarf galaxy population or, as we suggest, the 
turn up occurs at fainter magnitudes (ie. the 
dwarf to giant ratio, DGR, is smaller; we will de- 
fine the DGR as in Paper III as the number of 
galaxies with —16.5 > Mr > —19.5 compared 
to those with —19.5 > Mr). This later upturn 
is, in principle, verifiable with yet deeper pho- 
tometry, though background contamination un- 
certainties eventually dominate (see, e.g., Driver 
et al. 1997b = Paper II; Trentham 1997a). 

The clusters in question are not distinguished 
by their richness, but we can also check for 
morphological differences. As is well known, 
structural and (giant galaxy) population char- 
acteristics are well correlated, with, for exam- 
ple, dense regular clusters being of early Bautz- 
Morgan type (dominated by cD galaxies) and 
having the highest fractions of giant ellipticals. 
We can therefore choose to characterise the clus- 
ters by their central (giant) galaxy number densi- 
ties, for instance the number of galaxies brighter 
than Mr = —19.5 within the central 1 Mpc 2 
MI CM. (An alternative would be to use Dressler's 
(1980) measure of the average number of near 
neighbours, see Paper III). We find that the clus- 
ters with less prominent dwarf populations (lower 
DGRs ~ 1) are just those with the highest pro- 
jected galaxy densities (eg. A3888). The equiva- 
lent effect of earlier B-M type clusters having less 
dwarfs (A3888 is B-M Type I-II) was illustrated 
in Paper III (see also Lopez-Cruz et al. 1997). 
Earlier, Turner et al. (1993) had noted that the 
rich but low density cluster A3574, which is very 
spiral rich (Willmer et al. 1991), had a very 
high low surface brightness (LSB) dwarf to gi- 
ant ratio. This is now backed up by the obser- 
vations (see Paper III) of the clusters like A204 
which are dwarf rich (DGR ~ 3), have low cen- 
tral densities and late B-M types (A204 is B-M 
III). In addition, it appears that the Virgo Clus- 
ter, an archetypal loose (moderately rich) clus- 
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ter, has a very large dwarf population (Binggeli, 
Sandage & Tammann 1985), with many LSB 
dwarfs down to the sizes of Local Group dwarf 
spheroidals (Schwartzenberg, Phillipps & Parker 
1996; Phillipps et al. 1998). 

2.2. Population Gradients 

It was already suggested by the results on 
A2554 in Paper I, that the dwarf population was 
more spatially extended than that of the giants, 
ie. the dwarf to giant ratio increased outwards. 
This type of population gradient has been con- 
firmed by the results in Paper III. (We also con- 
sider this in more detail elsewhere, Smith et al., 
in preparation = Paper IV, using observations 
extending over larger areas). It is found, too, in 
Virgo (Phillipps et al. 1998), where the dwarf 
LSBG population has almost constant number 
density across the central areas while the giant 
density drops by a factor ~ 3. For Coma, the 
apparent discrepancy between, for example, the 
LF slopes of Bernstein et al. (1995), for the core, 
and Biviano et al. (1995), for a larger area, could 
similarly be attributed to an increase in the dwarf 
fraction in the outer parts (see also Thompson & 
Gregory 1993, Karachentsev et al. 1995, and sec- 
tion 3 below). 

2.3. The Field 

The opposite extreme of environment to the 
elliptical rich core of a dense cluster is, of course, 
the spiral rich 'field', which we can think of as 
made up of loose groups, and possibly the out- 
skirts of richer systems. Extending the above 
arguments we would therefore expect to see a 
very large dwarf population in the field. This 
is against the usual perception of the field LF 
at the faint end (e.g., Efstathiou, Ellis & Peter- 
son 1988; Loveday et al. 1992), but evidence has 
begun to accumulate that the field LF may in- 
deed turn up at the faint end, at a similar point 
to that seen in the clusters, around Mr ~ —19 
(e.g., Marzke, Huchra & Geller 1994; Zucca et 
al. 1997; see also Driver & Phillipps 1996). Very 



recently, considerations of the satellites of nearby 
spiral galaxies have suggested that the field dwarf 
LF may be just as steep as that in clusters (Mor- 
gan, Smith & Phillipps 1998), or, indeed, much 
steeper (Loveday 1997). A contributory factor to 
this revision of the field dwarf LF is clearly the 
inclusion of LSBGs missed from earlier surveys 
(Phillipps & Driver 1995; Ferguson & McGaugh 
1995), but it should be said that in the Local 
Group, where extremely low surface brightness 
dwarfs are, in principle, detectable from their re- 
solved stars, the LF appears quite flat. From the 
relationship derived in Paper III, a field LF with 
a faint end slope, say, a = —1.5 would have a 
DGR of 4 - 5. This would be nicely consistent 
with the extrapolation of the trend seen in the 
clusters (see below) down to low (volume) densi- 
ties. 

3. A Population Density Relation 

The obvious synthesis of the above results is 
to posit a relationship between the local galaxy 
density and the fraction of dwarfs (ie. the rela- 
tive amplitude of the dwarf LF). The inner, dens- 
est parts of rich clusters would have the smallest 
fraction of dwarfs, while loose clusters and the 
outer parts of regular clusters, where the density 
is lower, have high dwarf fractions. This is illus- 
trated in Figure 1, which is based on our homoge- 
neous rich cluster data from Papers I and III. It is 
particularly interesting to note the clear overlap 
region, where regions of low density on the out- 
skirts of dense clusters (open squares) have the 
same DGRs as the regions of the same density at 
the centres of looser clusters (solid squares) . The 
triangles show in slightly more detail the run of 
DGR with radius (hence density) across an indi- 
vidual cluster, A2554. 

The proposed relation of course mimics the 
well known morphology - density relation (Dressier 
1980), wherein the central parts of rich clus- 
ters have the highest early type galaxy fraction, 
this fraction then declining with decreasing lo- 
cal galaxy density. Putting the two relations to- 
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gether, it would also imply that dwarfs preferen- 
tially occur in the same environments as spirals. 
This would be in agreement with a weaker clus- 
tering of low luminosity systems in general (Love- 
day et al. 1995; see also Dominguez-Tenreiro et 
al. 1996) as well as for spirals compared to el- 
lipticals (Geller & Davies 1976). Thuan et al. 
(1991) have previously discussed the similar spa- 
tial distributions of dwarfs (in particular dwarf 
irregulars) and larger late type systems. 

In Figure 2 we have added to our data (squares) 
values derived from the work of other observers. 
A problem here is, of course, lack of homogene- 
ity due to different observed wavebands, different 
object detection techniques and so forth. Never- 
theless, we can explore whether literature results 
support our results. Firstly, several points are 
shown for various surveys of Coma (hexagons). 
These surveys (Thompson & Gregory 1993, Lobo 
et al. 1997, Seeker & Harris 1996 and Trentham 
1998) cover different areas and hence different 
mean projected densities. All these lie close to 
the relation defined by our original data, with the 
larger area surveys having higher DGRs. Points 
(filled triangles) representing the rich B-M type I 
clusters studied by Lopez-Cruz et al. (1997) fall 
at somewhat lower DGR than most of our clus- 
ters at similar densities. However we should note 
that these clusters were selected (from a larger 
unpublished sample) only if they had LFs well 
fitted by a single Schechter function. This obvi- 
ously precludes clusters with steep LF turn-ups 
and hence high DGR. The one comparison clus- 
ter they do show with a turn up (A1569 at DGR 
~ 4.2) clearly supports our overall trend. Al- 
though there is now considerable scatter (and the 
errors on some of the points are quite large), a 
weighted least squares fit to the trend gives log 
(DGR) = const. - (0.86 ± 0.22) log (giant den- 
sity), indicating a significant variation. Within 
the overall rich cluster sample, if we differentiate 
by the clusters' B-M types there is a suggestion 
that the type Is lie at lower DGR than the others, 
but again this may be biased by the Lopez-Cruz 



sample's selection criteria. Finally, a steep field 
LF (a ~ —1.5) would also give a point (filled 
pentagon) consistent with the trend seen in the 
clusters. 

On the other hand, Ferguson & Sandage (1991 
= FS), from a study of fairly poor groups and 
clusters, deduced a trend in the opposite direc- 
tion, with the early type dwarf to giant ratio in- 
creasing for denser clusters. However, this is not 
necessarily as contradictory to the present result 
as it might initially appear. For instance, FS 
select their dwarfs morphologically, not by lumi- 
nosity (morphologically classified dwarfs and gi- 
ants significantly overlap in luminosity) and they 
also concentrate solely on early type dwarfs. If, 
as we might expect, low density regions have sig- 
nificant numbers of late type dwarfs (irregulars), 
then the FS definition of DGR may give a lower 
value than ours for these regions. Furthermore 
FS calculate their projected densities from all de- 
tected galaxies, down to very faint dwarfs. Re- 
gions with high DGR will therefore be forced to 
much higher densities than we would calculate 
for giants only. These two effects may go much 
of the way to reconciling our respective results. 
This is illustrated by the open triangles in figure 
2 which are an attempt to place the FS points on 
our system; magnitudes have been adjusted ap- 
proximately for the different wavebands, DGRs 
have been estimated from the LFs and the clus- 
ter central densities (from Ferguson & Sandage 
1990) have been scaled down by the fraction of 
their overall galaxy counts which are giants (by 
our luminosity definition). Given the uncertain- 
ties in the translation, most of the FS91 points 
now lie reassuringly close to our overall distribu- 
tion. 

Nevertheless, there are two exceptions, the FS 
points of lowest density (the Leo and Doradus 
groups) which also have low DGR (and lie close 
to our main 'outlier', the point for the outer re- 
gion of A22). The Local Group (shown by the 
star) would also be in this regime, at low density 
and DGR = 2, as would the 'conventional' field 
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with a ~ —1.1 and DGR ~ 1.5 (open pentagon). 
This may suggest that at very low density the 
trend is reversed (ie. is in the direction seen 
by FS), or that the cosmic (and/or statistical) 
scatter becomes large. More data in the very 
low density regime is probably required before 
we can make a definitive statement on a possible 
reversal of the slope of the DGR versus density 
relation. In particular, the scatter in the derived 
faint end of the field LF between different sur- 
veys (see, e.g., the recent discussion in Metcalfe 
et al. 1998) precludes using this to tie down the 
low density end of the plot. 

4. Discussion and Summary 

As with the corresponding morphology den- 
sity relation for giant galaxies, the cause of our 
population - density relation could be either 'na- 
ture' or 'nurture', ie. initial conditions or evolu- 
tion. Some clues may be provided by the most 
recent semi-analytic models of galaxy formation, 
which have been able to account in a general way 
for the excess of (giant) early type galaxies in 
dense environments (e.g., Baugh, Cole & Frenk 
1996). 

The steep faint end slope of the LF appears 
to be a generic result of hierarchical clustering 
models [] (e.g., White & Frenk 1991; Frenk et al. 
1996; Kauffmann, Nusser & Steinmetz 1997), so 
is naturally accounted for in the current genera- 
tion of models. The general hierarchical forma- 
tion picture envisages (mainly baryonic) galaxies 
forming at the cores of dark matter halos. The 
halos themselves merge according to the general 
Press-Schechter (1974) prescription to generate 
the present day halo mass function. However 
the galaxies can retain their individual identi- 
ties within the growing dark halos, because of 
their much longer merging time scales. The ac- 
cretion of small halos by a large one then results 
in the main galaxy (or cluster of galaxies, for very 

1 And was considered a problem until observational evi- 
dence for steep LFs increased! 



large mass halos) acquiring a number of smaller 
satellites (or the cluster gaining additional, less 
tightly bound, members). 

Kauffmann et al. (1997) have presented a de- 
tailed study of the distribution of the luminosi- 
ties of galaxies expected to be associated with a 
single halo of given mass. The LFs are some- 
what disjoint owing to the specific halo masses 
modelled; especially for the low mass halo there 
is a preferred luminosity for the central galaxy 
plus a tail to lower luminosities. For a realis- 
tic mix of halo masses, these would no doubt be 
smoothed to look more like conventionally ob- 
served LFs. Nevertheless, we can still easily com- 
pare the numbers of dwarf galaxies per unit giant 
galaxy luminosity (rather than the amplitude of 
the giants' LF) between halos of different mass. 

The Kauffmann et al. models mimic a "Milky 
Way system" (halo mass 5 x 10 12 Mq), a size- 
able group (halo mass 5 x 10 13 Mq) and a clus- 
ter mass halo (10 15 Mq). Using their figure 2 
(which also emphasises the identical faint end 
slopes predicted for all the different environ- 
ments), we choose to quantify the number of 
dwarfs by iV_is, the number of dwarfs per system 
in the Mb = — 18 bin. Because of the very similar 
slopes, the choice of bin or range of bins does not 
affect our conclusions, so this is the equivalent of 
the total number of dwarfs used in Figure 1. To 
quantify the giant population we choose the to- 
tal light of galaxies of Mb = —20 or brighter, in 
units of galaxies (taking M B = —21), which 
we call N_2i- Using this definition, rather than 
the actual number of galaxies brighter than some 
value (as in our observational data) allows for the 
discretization of the LFs for small halos. The re- 
sults are summarized in Table 1. The ratio of 
these two values iV_i8 and N_2i then quantifies 
the relative dwarf galaxy populations. Roughly 
speaking, for smooth LFs with a shape similar 
to that observed, we should multiply these val- 
ues by about 5, giving a range from about 1 to 
3, to compare with our observational DGRs. 

We see that the Milky Way and small group 
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Halo Mass (M ) 


iV_ 18 


N-21 


iV_ 18 /iV_ 21 


5 x 10 i2 


0.2 


0.34 


0.6 


5 x 10 13 


2.2 


3.8 


0.6 


1 x 10 15 


40 


190 


0.2 



Table 1: Dwarf numbers as a function of halo mass. 

halos have similar numbers of dwarf galaxies per 
unit giant galaxy light, whereas the dense clus- 
ter environment has a much smaller number of 
dwarfs for a given total giant galaxy luminosity. 
Thus the predictions of the hierarchical models 
(which depend, of course, on the merger history 
of the galaxies) are in general agreement with our 
empirical results if we identify loose clusters and 
the outskirts of rich clusters with a population of 
(infalling?) groups (cf. Abraham et al. 1996), 
whereas the central dense regions of the clus- 
ters originate from already massive dark halos. 
By inputting realistic star formation laws etc., 
Kauffmann et al. can further identify the galax- 
ies in the most massive halos with old elliptical 
galaxies, and those in low mass halos with galax- 
ies with continued star formation. This would 
imply the likelihood that our dwarfs in low den- 
sity regions may still be star forming, or at least 
have had star formation in the relatively recent 
past (cf. Phillipps & Driver 1995 and references 
therein). Note, too, that these galaxy formation 
models would also indicate that the usual (giant) 
morphology - density relation and our (dwarf) 
population - density relation arise in basically the 
same way. Finally, we can see that if these mod- 
els are reasonably believable, then we need not 
expect the field to be even richer in dwarfs than 
loose clusters; the dwarf to giant ratio seems to 
level off at the densities reached in fairly large 
groups. 

To summarise, then, we suggest that the cur- 
rent data on the relative numbers of dwarf galax- 
ies in different clusters and groups can be un- 
derstood in terms of a general dwarf population 
versus local galaxy density relation, similar to 



the well known morphology - density relation 
for giants. Low density environments are the 
preferred habitat of low luminosity galaxies; in 
dense regions they occur in similar numbers to 
giants, but at low densities dwarfs dominate nu- 
merically by a large factor. This fits in with the 
general idea that low luminosity galaxies are less 
clustered than high luminosity ones (particularly 
giant ellipticals). Plausible theoretical justifica- 
tions for the population - density relation can be 
found within the context of current semi-analytic 
models of hierarchical structure formation. 
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Fig. 1. — Variation of the dwarf-to-giant ratio (DGR), as denned in the text, with projected density of 
cluster giants (per square Mpc). Solid boxes represent the central 1 square Mpc regions of the clusters, 
the open boxes the outer regions (data from Paper III). The triangles show the variation over a wider 
range of radii for Abell 2554 (data from Paper 1). Note that typical error bars (due to the combination 
of poisson errors and background subtraction errors) are 10% in density and 20% in DGR for the denser 
regions, rising to 30% in density and 50% in DGR at the lowest densities (and hence object numbers). 
The outlier at low density and low DGR (the outskirts of A22) has a very large error in DGR (~ 100%). 
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Fig. 2. — As figure 1, but including data from other observers. Squares are our data repeated from figure 
1, hexagons are for various Coma surveys detailed in the text, filled triangles are from Lopez-Cruz's 
sample and open triangles are for Ferguson and Sandage's poor clusters and groups. The open pentagon 
at low density represents a conventional 'flat' field LF, the filled pentagon a possible steep (a ~ —1.5) 
field LF and the Local Group is represented by the star at DGR = 2. 
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